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Abstract 

The square integrable basis set representation of the resolvent of the asymptotic three-body 
Coulomb wave operator in parabolic coordinates is obtained. The resulting six-dimensional Green's 
function matrix is expressed as a convolution integral over separation constants. 
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I. INTRODUCTION 



It is well known that the Schrodinger equation for a three-body Coulomb system is 
asymptotically separable in terms of parabolic coordinates [l| 



£j = r ls + k /s • r ls , rjj = r is - k ls ■ r /s 



(1) 



where r/ s and are the relative coordinate and momentum vectors between the particles 
/ and s. Here j, I, s is a cyclic permutation of 1, 2, 3. The long-ranged six- dimensional 
operator, which provides a three-body continuum wave function, 3C-function 

as 

, that 

satisfies exact asymptotic boundary conditions for Coulomb systems in the region where the 
distances between the particles are large, reads as the sum of two-dimensional operators 

3 

i 

h 3 + h m + %hsU 



y L . 



1 : 



(2) 



where t ls = fi ls 
by 



minis 
mi+m, 



. Here the one-dimensional operators and h Vj are defined 



h 



(3) 



In the previous paper 4| we introduce the six- dimensional operator f) which is obtained by 
multiplying (j2J) on the left by H j=1 +Vj)'- 

h = ^13 (6 + V2) A*12 (& + 773) *)l + /i 2 3 (6 + ^1 ) ^12 (6 + *7s) 62 + /i23 (6 + »7l ) ^13 (6 + ffe) (>3 , (4) 

where 

§3 = hj + Kj + Ihstla- (5) 

The resolvent of the operator (Tj0) can be used in the corresponding Lippmann-Schwinger 
equation for the three-body Coulomb wave function. 

It has been suggested [4( to treat the operator (j4]) within the context of L 2 parabolic 
Sturmian basis set [5| 

3 

(6) 



3=1 

^ntrnjiXj, Vj) = '•'» i (s / ) >■;„,{ I/ j). 



2be~ bx L n (2bx), 



(7) 
(8) 



where b is the scaling parameter. In particular, a matrix representation Gj of the resolvent 
for the two-dimensional operator 

h, • fnJ'jiC, ■ ni) (9) 

has been obtained which is formally the matrix inverse to the infinite matrix [hj + fii s CjQ,j] 
of the operator (J9]): 

[hj + MsCjQj] Gj(t ls ; fiisCj) = Ij. (10) 

Here 

hj = <g> I Vj + I 5j ® h v . + 2ki s ti s \j (11) 

is the matrix of the operator ()j ([5]) in the basis ([7]), 1^., I % and Ij = 1^. <8> I Vj are the unit 
matrices. Qj = Q^. <g> 1^ . + 1^. ® Q^. , where Q^. and are the matrices of £j and in 
basis (jHJ), respectively. 



In the previous paper J4] it has been suggested that the matrices Gj of the two- 
dimensional Green's functions be used to construct the matrix (5 which is inverse to the 
operator (j3J) matrix 

f) = ^la/^uhi <8> Q 2 ® Q 3 + M23AH2Q1 ® h 2 <g> Q 3 + /i 23 /ii3Qi ® Q2 ® h 3 . (12) 

Namely, we proposed to express the six- dimensional Green's function matrix (3 as the con- 
volution integral 



= K / / dC lC ?C 2 G!(t 23 ; #bA) ® G 2 (t 13 ; // 13 C 2 ) ® G 3 (t 12 ; -^12(^1 + C 2 )), (13) 




where K is a normalizing factor. Thus, our problem is now to determine the pathes of 
integration over the separation constants C\ and C 2 in ffl3|) and to find the corresponding 
normalizing factor K such that the condition 

1)0 = Ii ® I 2 ® I 3 (14) 



3 



holds. For this purpose consider the product f) (3. From the relation (jTUl) we obtain 
h_0 = N{J7 dC 1 dC 2 [I 1 -fi 23 C 1 Q 1 G 1 (t 23 ; /i 23 d)]®^ 13 Q 2 G 2 (* 13 ; /i 13 C 2 ) 

®/Ui2Q 3 G 3 (ti 2 ; -Aii 2 (Ci + C 2 )) 
/ / dCidC 2 /i23QiGi(t23; [I 2 -^i 3 C 2 Q 2 G 2 (t 13 ; // 13 C 2 )] 

(15) 

®/Ui 2 Q 3 G 3 (ti 2 ; -/ii 2 (d + C 2 )) 
/ / rfCirfC 2 // 2 3QiGi(t 23 ; /X23C1) ® /Ui3Q2G 2 (ti 3 ; /ii 3 C 2 ) 
® [I 3 + ^i 2 (Ci + C 2 )Q 3 G 3 (t 12 ; -mi 2 (Ci + C 2 ))]} , 

and hence 

&0 = N{/ / dC lC ?C 2 [Ii ® ^i 3 Q 2 G 2 (t 13 ; /ii3C 2 )®/ii 2 Q 3 G 3 (t 12 ; -yL 12 {C x + C 2 )) 

+yU 23 QiG!(t 23 ; /123d) ® I 2 ® / u 12 Q 3 G 3 (t 12 ; -^(Ci + C 2 ))] (16) 

+ [//23Q1 / rfCiGi(t 23 ; imCx)] ® [/iiaQa / dC 2 G 2 (ti 3 ; //i 3 C 2 )] ® I 3 } . 
As a first step towards our goal we consider the integrals 

J dCjGfiu; fUsC,) (17) 

inside the figure brackets on the right-hand side of ( [161) . 

In Sec. II completeness of the eigenfunctions of the two-dimensional operator is con- 
sidered. In particular, an integral representation of the matrix A which is inverse to the 
infinite matrix Q of the operator (£ + rj) in the basis ([7]) is obtained. In Sec. Ill it is demon- 
strated that the integral (fTTI) taken along an appropriate contour is proportional to the 
matrix A obtained in the previous section. Finally, Sec. IV presents a convolution integral 
representation of the six-dimensional Coulomb Green' function matrix. 

II. COMPLETENESS RELATIONS 
A. The continuous spectrum 



Of particular interest are the regular solutions 

f(l, t , £, V) = w(t> r, i>(7, r, 77) 
4 



(18) 



of the system 



hz + 2kt + vC£ w( 7 , r, = 0, 
h v + 2k(t -t) + fiCri] u(7, r, = 0. 



(19) 
(20) 



Obviously, the regular solutions u and v are proportional to confluent hypergeometric func- 
tions 



and 



u( 7 , r, = e*^*» ^ Q + <t, 1, -ryf) 
v( 7 , r, rj) = e-fr-Qv lFl {\ + i( T - r ), 1, 7777) 



(21) 



(22) 



where 



t + 



k 

TO = -to- 

7 



(23) 



Y k 
2 _ T' T ~ 7 

It should be noted that since the representation of the two- and six-dimensional Coulomb 
Green's functions matrix elements involves an integration over r from — oo to oo, in the 
subsequent discussion we assume that r is real. With this assumption it is readily seen that 
the solutions (j2ip and fT22|) coincide, except for normalization and the phase factors e~5 fc £ 
and e^ kv , with parabolic Coulomb Sturmians treated in Ref. Q]. In this case 7 plays the 

2 

role of the momentum and £ = \ is the energy. From this we conclude that for 7 > the 
solutions (12ip and (1221) correspond to the continuous spectrum of £. 

It is readily verified that the solutions u(j, r, £) and v (7, r, 77) are expressed by basis set 
(JHJ) expansions 



w(7> r > 



2V26 



u( 7 , r, 77) 
where 



26- 1(7 - fc) 
26-i( 7 -fc) V.2& + ^(7 + fc) 

2^ 



00 



£>p n (r; C) V>n(£), 

n=0 



26 - ^(7 + fc) \2b + i(-f- k) 
2b + z( 7 - fc) 



26-i( 7 + ife)y('»-)+> « 

1 2^ A >n(r -r; C) V»nfa), 



n=0 







A 



2b - 1(7 + fc) 



2&-i( 7 -Jfe)' 26 + i( 7 + fc)' 

The expansion (1241 and (1251) coefficients contain the polynomials 



Pn( T ! C) 



-l) n T (n+ I -ir) 
r (I - ir) 



1 • 1 ■ , 

-n, - + -n + - + 2r; ( 



(24) 



(25) 



(26) 



(27) 
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The basis set © representation of the equation (fl9l) is the three-term recursion relation 

1 

a n Vn-i + b n y n + d n+ i y n+l = 0, n > 1 (28) 

where 

b n = (b+i§ + ik)+ 2(b + t§)n + 2kt, 

(29) 

a n = ( b - *w ~ ik ) n , d n = (b - + ik)n. 

The functions 

s n (t; fiC) = 6 n p n (r; () (30) 

are the "regular" solutions of ( |28i) with the initial conditions: so = 1, s_i = 0. The 
polynomials p n (T27j) of degree n in r are orthogonal with respect to the weight function {4] 

^0 = r({ '^ (i + l,1 (^', (3D 
where it is considered that |arg(— £)| < 7r. The corresponding orthogonality relation reads 



00 

^ (^(~) / dr P( T > OPn(r; C)Pm(r; C) = 5 nm . (32) 



B. The discrete spectrum 



For t < the eigenf unctions of (TT9~j) corresponding to the discrete spectrum: £^> = 
7f = ^ = -^a, £ = 1, 2, . . . , 00, are [8j 

fe, m {£, V) = ue,m(O v e,£-m-i(v), m = 0, 1, 1, (33) 

where 

w/,m(0 = e-5 fc «e~^ iFiC-m, 1; = e^e - ^ L m («/0 (34) 



and 

^.mfa) = e5*"e~T iFi(-m, 1; /c^) = e^e"^ L m {n e r)). t:V>) 
The solutions m meet the orthogonality relation 



00 00 

,3 



J{^ + V)d^dr]f etm (^,r])[f l ^ m ,(^,ri)]* = 5 e j l 5 rntm/ . (36) 
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It is readily verified that the expansions of U£ >m (£) an d v i,m(v) i n a basis function 
series are 



and 



n=0 



n=0 



where the coefficients are given by 



S, 



(l, m) 



2V2b 



\n (m+l) n (2b-K e -ik) n {2b~Kt+ik) r ' 
' n! (2b+Ki+ik) n + m + 1 



x 2 F x ( -n, -m, -n - m; ; 2fe _ K ;j 2+A;2 



(37) 



(38) 



(39) 



One-dimensional completeness relations 



The eigensolution u(j, r, £) of (|T9l corresponding to the continuous spectrum (7 > 0) 
for large £ behaves as 



e-»« r(i + zr) Vt£ V 2 



where a = argT (~ + it). Therefore (see e. g. [9j) 



e 2 fcf 2z . {i£ vr 

sin r m(74J H h a 



4 



(40) 



2 
4 

and for r > 



r.- + /r 



e^ r y e^«( 7 , r, o [«(Y, r, 0]* = ^(7 - 7), 





(41) 



£ / 7^7 



47T 



•tt( 7 , r, K 7 , r, = £')■ 



(42) 



On the other hand, if the functions ^(7, r, £) are regarded as charge Sturmians i. 
e. the parameter r is considered as eigenvalue of the problem, whereas the momentum 7 
remains constant, the corresponding orthogonality and completeness relations are given by 

a 



7 



' " I d^u^, t, ^{ui'j, t', = 2tt5(t - t') 




and 



7 / dr 



2tt 



■e— u( 7 , r, 0K 7 , r, Of = O- 



(43) 



(44) 



7 



Taking matrix elements of the completeness relation (1441) we find 



oo 



' dr i—^- — (-C) Pn (r; C) [Pn (r; C)] = ^ 



v/[46 2 + (7 + k) 2 } [4b 2 + (7 - A;) 2 ] 7 2vr 

—00 

(45) 

It may be noted that is closely related to (f3"2"l) . To see this, let 7 > 0. Then, it follows 
from the definitions (1261) that 

^(-0* - (46) 
C V[46 2 + (7 + k) 2 } [46 2 + (7 - k) 2 } 

Further, the regular solution 8 n p n (r; Q of the three-term recursion relation (1281) is an even 
function of 7, since the coefficients a n , b n and d n depend on 7 only through fiC = \{k 2 — 7 2 ). 
Thus, replacing 7 by —7, and hence r by — r (0 — > A and A — > 0) and ( by in Eq. ( 1271) 
gives 



9 n Pn (r; C) = A» ( J ! )Wr ^ ( t|^ r) 2 F 1 (-n,l-.r; -n + ± - ir; C" 1 )- (47) 
Comparing Eqs. (f2~Tj) and fj471) then yields the relation 

#>n(r; C) = A n b„(r; C)]* , (48) 

and hence 

\Pn(r; Or = r n Pn(r; 0- (49) 



From the argument above, we conclude that for 7 > the orthogonality relation (1321) reduces 
to 



The two-dimensional completeness relation 

It follows from the relations (141D and (14*31) and analogous relations for v (7, r, £) that the 
two-dimensional orthogonality relation for 7 > is given by 



y 2 l r Q+^)l I r (l+^Q-^))| 

4 2?r 2tt 



2 00 00 



e r rvj^vi ru T -vi / /(e + 77) dedf7 {/(t, r, C, v) 



(50) 
x[/(Y, r', r ? )r} = 5( 7 -7 / )^-r'). 
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In turn, in the case to > (where there are no bound states) it would appear reasonable 
that the two-dimensional completeness relation would be given by 

i r Q+^)ri r Q+^o-^))i 2 



( oo oo 

(f + r]) I a J d>y>y 2 e- 7TT0 J dr l - 

l -oo 



2vr 2vr 

(51) 

x/(7, r, £, V) lf(l, t, C v')Y \ = 6(t- O 6(r) - rf). 



The integration over r in (15T1) is performed on the assumption that r is independent of 7. 
To test this hypothesis and determine the normalizing factor a, we carried out the following 
numerical experiments. First with some parameters to > 0, k > and b we calculate the 
matrix elements A nijn2;mijm2 for the expression in the figure braces on the right-hand side 
of ([EED in the basis (0): 



A = nfd^V 64 b 2 7 2 (-Cr" 0ni-mi\ma-n 2 

/1 "l,"2;mi,m2 "J" / f4& 2 + (^+M 2 lf4& 2 +^-fc) 2 l 



X 
—00 



00 2° 2 (52) 

1 dr |r( ^ T)|a|r( ^" T))|a p ni (r; C)Pn 2 (r - r; C) K (r; C)p m2 (r - r; QT ■ 



It should be noted that the value of (— (Y T ° in this formula is determined by the condition 
|arg(— C)| < 71. Then the resulting matrix A is multiplied by the matrix 

Q = Q 5 <g> l r , + I c ® (53) 

of the operator (£ + rj). Finally, using the condition 

QA = I^I T (54) 

we have obtained that a = \. Notice that the infinite symmetric matrices and Q,, are 
tridiagonal 4|, therefore the equations on the left hand-side of the linear system (1541) each 
contain no more than six terms. 

For for to < the completeness relation ( 15T|) transforms into 

(Z + V) {|7rf 77 2 e— J dr ^ T){2 |r(l+ tT T))|2 /( 7> r, £, 17) [/( 7 , r, 

I ° 3 « ~°° 1 (55) 

£=1 m=0 J 

In this case the matrix A with elements 

A f^- 32fe 2 7 2 (-C)' T o nn 1 -m 1 \m 2 -n 2 

^n 1 ,n 2 ;m 1 ,m 2 J a J [ 46 2 + ( 7+fc )2] [ 4fe 2 + ( 7 _ fc )2] P ^ 

x 7 rfr |r(l 2 + ; r)|2|ra+ t°" T))|2 ^ CK 2 (to - r; C) [p Wl (r; C)Pm 2 fa - r; C)]* (56) 



-00 

00 , t-\ 



k£ c(^' m ) c(^' ^ — m— 1) 

2^ 2I 2^ ni m2 

£=1 m=0 



is also inverse to the matrix Q (|53|) . The expression (1561) can be rewritten, in view of (146 
and (TJ9]), as 



A 



ni, n2; mi , m2 



I v 7 -00 
XPni(T] C)Pn 2 (r -T; C)p mi (r; C)p m2 (r -T; () 

00 o £-1 



•>{l,m) q(£,e-m-l) 

=1 m=0 



o(£,£-m-l) q(e,m) 



(57) 



III. CONTOUR INTEGRALS 



g n (r, C) - {-) r(n+f+ir) 2*1 U + 



(fSrV; = (-)' 

In particular, the functions 

7(±) 



»!rU-*rl 



2 F X (I — it, n + 1; n + § - ir; () . 



w (±) (7, 7", 



2iyJL 



2b+i( 7 -fc) y r+ 2 



26-i(7-fc) ^ 26-j(7+fc) y 



V 2 I n=0 



tend to w^(j, r, £) as £ — > 00. 

The matrix elements of the resolvent of 



9n 

(t; (mC) 



+ 2kt + fx C£ 
C — 1\ 9" 



(58) 



Notice that expressing the resolvent of the one- dimensional operator + 2kt + \i C£ 
requires two linearly independent solutions of f[T9l . Irregular solutions of (ITU1) are expressed 
in terms the confluent hypergeometric functional 

W W( 7 , r, = ei^-^[/ Q ±<t, 1; pye). 

The corresponding solutions of the three-term recursion relation (128]) are 

ci+)(t ; M C) = r +1 (r; C) , 
d~\t; yU C) = A"+ 1 g(- ) (r; C), 

where 

,C+)^. ^_^n^G+-) ^ n , jTjn + 1 . n+ l + lT . C -i) ) 



(59) 



(60) 



(61) 



can be written in the form 



4] 



and 



27V c ; c 

C - 1 \ 9 n 



P n< (r; OdtV; 



27V c ; c 



Pn< (r; C) C >+1 qL\r; Q, 



(62) 



(63) 
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where n> and n < are the greater and lesser of n and m. Notice that c« + ' ) ( c «~' ) ) are defined 
in the upper (lower) half of the complex 7-plane where |C|>1(|£|<1). To analytically 
continue onto the lower half of the 7-plane the relation (4] 



c«(£; /iC) = c(-)(t; pC) + 2mp(r ] () 9 n+1 p n (r; 



(64) 



can be used. 

In {4] we obtained the basis set ([7]) representation of the resolvent for the two-dimensional 



operator 



+ 2kt + fiCl; + h v + 2/c(to — t) + /iCr^ . In particular, the matrix elements 



of the two-dimensional Green's function can be expressed as the convolution integral 

G nitn2] mi , m2 (t ; fiC) — i (—^j A 2 / dr p(r — r; Q 

—00 

x&ii*; ^C) p n2 (T - t; ()p m2 (T -T; C). 



(65) 



Notice that in this case only the regular solutions of (|20|) discrete analogues A~ n p n (r — r; 
are used. 

Let us consider the integral 



(66) 



i/C, 



T 



and 



Notice that by replacing 7 — > —7 (and hence 6 — ► A, A — > 0, £ 
To * To) in Eq. (1651) Gn!,^; mi,m2 (t ; yuC) is transformed to Gni,n 3 ; m 1: m 2 {to] pC). Thus, 
for the integral X\ we obtain 



00 



(67) 



Inserting Eqs. (16"5"j) . (1621) and (1651) into Eq. (157)) . we find, in view of (I54"j) . that Xi coincides 
with the integral on the right hand-side of (l5"7j) . 
Now, we consider the integral 



— / d£G« f t ; 
2vri / 1 2 



(68) 



taken along the contour in the complex £-plane shown in Fig. 1. The contour C passes in a 
negative direction (clockwise) round all the points £^> = — ^ (filled circles in Fig. 1 which 
accumulate at the origin) and the cut along the right half of the real axis and is closed at 



11 



infinity (see e. g. [id]). The corresponding matrix element of the integral along the two sides 
of the cut is equal to X\ f l66|) (this is circumstantial evidence that the normalizing factor a 
in the completeness relation (1511) is equal to |). On the other hand, the integration along a 
contour enclosing £^> reduces to (—1) times the double sum of the residues of the integrand 
at the points t^> — i (m + |), m — 0, 1, ... and £^ = ; £ — 1, 2, ... , which are 

the poles of the gamma functions F (| + ir) and r (| + i(r — r)) = T ym + 1 + —j ^= j , 
respectively. It is readily shown that the matrix element of this part of the integral ( IBS]) 



coincides with the double sum in floTj) . Thus, the integral fl68|) is equal to the matrix A. 
The contour C can be deformed, for instance, into a straight line parallel to the real axis. 
The resulting path C\ shown in Fig. 2 runs above the cut and the bound-state poles of 
G^" 1 ") (to; % — £ ). Then, to make the integral amenable to numerical integration, rotate C\ 

n 

about some point £0 on the right half of the real axis through a negative angle [10J; see 
the contour C 2 in Fig. 2. The part of C2 on the unphysical sheet is depicted by the dashed 
line. Thus, we obtain the following representation of the matrix A (which is inverse to the 
matrix Q f l53|) ): 

A = ^/ d£G<+> ( t » i T- £ )= Q " (69) 

IV. SIX-DIMENSIONAL GREEN'S FUNCTION MATRIX 

Using the relation (|69l) we can rewrite the expression (fT3|) for the six-dimensional Coulomb 
Green' function matrix as the contour integral 

1 • - — — 23 c Ln * .13 c \ o /~i / + . 12 



- = l^m J d£l J d£2 Gl V 23; 2 ~ 7 ® ° 2 \ 135 2 ~ f2 J ® ° 3 V 12; 2 " 83 

c 2 c 2 

(70) 

fc 2 

where £j = ni s Cj and C1+C2+C3 = 0. This also allows us to determine the normalizing 
factor K. Indeed, it follows from ( 1691) that the third term inside the figure brackets on the 
right-hand side of f|T6l) is proportional to the unit matrix: 



Qx / dS x Gi (t 23 ; % - £t) ® Q 2 / d£ 2 G 2 (* 13 ; f - £ 2 
c 2 v J c 2 v 



Is 

(71) 



(27Ti) 2 l! ® I 2 ®L 
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Consider the first two terms in the figure braces in (fTBj) . For the energy 8 3 = -^+^12(^1+^2) 
we have 



1^12 ( k 



2 /i 23 V 2 



/"12 / 



A*i3 V 2 



"13 



On the other hand, on the contours C2 the energy variables £7, j — 1, 2 are given by 



(72) 



(73) 



where (p < 0, £oj is an arbitrary positive parameter, i?^ is real and runs from —00 to 00. 
Hence, for the energy £ 3 dZ2D we obtain 



kf 2 1^12 ( k 



2 /i23 V 2 



A*i3 V 2 



"13 



--02 



Thus, £3 can be expressed in the form 



£3 = £03 + E 3 e 



up 



m^-^Eh)**. (74) 

^23 1^13 



(75) 



where £03 an d E 3 denote the term in the square braces and the real factor in front of the 
exponent in (IT4"j) . respectively. Since £ 03 should be positive, the positive parameters £ i and 
£02 have to satisfy the constraint 



^12 c , ^12 c ^ ^12 , Ml2 ^23 , Ml2 &13 

C01 H 002 < — -\ — H — • 

/^23 A*13 2 /^23 2 yUl3 2 



Now, we consider the integral 



J 2 = / d£ 1 G 3 I t 12 ; -^-£3 



c 2 



(76) 



(77) 



With fixed £ 2 , in view (CU), ([75]) and 0, we see that 

1 2 = e* f dE 1 G 3 (t 12 ; k -f - £ 3 ) = e* / dE 3 G 3 Ui 2 ; k -f - £ 

—00 —00 

= "SJ ^3G 3 (f la ; f - £3) = ~2m g Q 3 \ 

Similarly, we obtain 

| d£ a G s ^ - £3) = -2*1 ^ Q-i. 
c 2 

Inserting (|75jl . (j79"j) and (|7Ijl into (USD then yields 

f> = K4tt 2 Ii ®I 2 ®I 3 - 



(78) 



(79) 



Thus, from (IHUD we conclude that 



4vr 2 



(80) 



(81) 
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V. CONCLUSION 



The Sturmian basis-set representation of the resolvent for the asymptotic three-body 
Coulomb wave operator is obtained, which can be used in the discrete analog of the 
Lippmann-Schwinger equation for the three-body continuum wave function. The six- 
dimensional Green's function matrix is expressed as a convolution integral over separation 
constants. The integrand of this contour integral involves Green's function matrices corre- 
sponding to the two-dimensional operators which are constituents of the full six-dimensional 
wave operator. The completeness relation of the eigenfunctions of these two-dimensional op- 
erators is used to define the appropriate pathes of integration of the convolution integral. 
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Figure 2: The path of integration C2 in (|69|) . 
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